Dynamics of a trapped Brownian particle in shear flows 
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The Brownian motion of a particle in a harmonic potential, which is simultaneously exposed ei- 
ther to a linear shear flow or to a plane Poiseuille flow is investigated. In the shear plane of both 
flows the probability distribution of the particle becomes anisotropic and the dynamics is changed 
in a characteristic manner compared to a trapped particle in a quiescent fluid. The particle distri- 
bution takes either an elliptical or a parachute shape or a superposition of both depending on the 
mean particle position in the shear plane. Simultaneously, shear-induced cross-correlations between 
particle fluctuations along orthogonal directions in the shear plane are found. They are asymmetric 
in time. In Poiseuille flow thermal particle fluctuations perpendicular to the flow direction in the 
shear plane induce a shift of the particle's mean position away from the potential minimum. Two 
complementary methods are suggested to measure shear-induced cross-correlations between particle 
fluctuations along orthogonal directions. 

PACS numbers: 87.15.Ya,05.40.-a,83.50.Ax 



I. INTRODUCTION 

The Brownian motion of particles in a fluid is of central 
importance in chemical and biological physics as well as 
in material science and engineering |l|-l4(. Despite the 
long history of Brownian motion, especially in quiescent 
fluids, our understanding of thermally induced particle 
dynamics in flows is still far from complete. 

Moreover, neutral colloidal particles moving relatively 
to each other interact via the fluid and these hydrody- 
namic interactions can cause a complex collective behav- 
ior 0, [H, In shear flows little is known about the 
dynamics of Brownian particles and the hydrodynamic 
interaction effects in spite of their fundamental relevance 
and importance in microfluidic applications. The Taylor 
dispersion [?} and fluid mixing issues 0, H[ are well known 
examples where fluctuations of particles and their hydro- 
dynamic interaction effects in simple shear and Poiseuille 
flow play an important role. The interplay of shear gra- 
dients and thermal motion of polymers leads, even at low 
values of the Reynolds number, to the rich dynamics of 
polymers [|[, the so-called molecular individualism [Io| . 
Polymers tumbling in a shear flow cause elastic turbu- 
lence even in diluted polymer solutions [11] and spectac- 
ular mixing properties [8| in microchannels. 

In shear flows it is the contribution (u • V)u of the flow 
field u(r) to the Navier-Stokes equation which causes in- 
teresting transient phenomena near the onset of turbu- 
lence |12| , as well as amplifications of velocity fluctua- 
tions and their cross-correlations along and perpendicu- 
lar to straight streamlines [l3|, Q • A cross-correlation is 
also expected between orthogonal particle- fluctuations in 
the shear plane, because random jumps of a particle be- 
tween neighboring streamlines of different velocity lead to 
a change in the particle's velocity and displacement along 
the streamlines, similar as via fluctuations. Nevertheless, 
there was no direct observation of these cross-correlations 
until recently [HI . Here we present the theoretical back- 



ground for their determination. 

The stochastic dynamics of free single spherical parti- 
cles in linear shear flows have been studied in terms of the 
hydrodynamic fluctuation theory [16], and in combined 
Langevin and Fokker-Planck approaches, by taking iner- 
tia into account [l7l-[l9| . Even effects of non-equilibrium 
thermodynamics were included in Ref. [2(j. Experi- 
ments for detecting shear-induced cross-correlations be- 
tween perpendicular random displacements of free parti- 
cles were described in Ref. [2l|, [22j . Shear induced cross- 
correlations between perpendicular fluid- velocity fluctu- 
ations and perpendicular fluctuations of particles are ex- 
pected to be strongly asymmetric in time [l3|, [HI, |23| . 
In dynamic light-scattering experiments certain aspects 
of this issue were observed indirectly [24|, however, a di- 
rect measurement and characterization of related particle 
fluctuations in shear flows remained an open question. 

Direct observations of particle fluctuations at the 
mesoscale became possible only quite recently by using 
optical tweezers. This rather young technique is a pow- 
erful experimental method for investigating the motion 
of a small number of particles [25l-[27| , which contributes 
substantially to our understanding of the dynamics of 
particles and to a number of innovative applications. 
These include the inspiring studies on single polymers 
[28l-[35j, the detection of anti-correlations between hy- 
drodynamically interacting Brownian particles by Femto- 
Newton measurements [36j , the propagation of hydrody- 
namic interactions [37j , wall effects on Brownian motion 
38 , [39| , short-time inertial response of viscoelastic fluids 
|40| , two-point microrheology [I]] , anomalous vibrational 
dispersion [42| , particle sorting techniques (43l-[46| , and a 
number of further investigations in microfluidics. The 
laser tweezer technique has also been applied to deter- 
mine the force elongation relation of biopolymers 1471 or 
the effective pair potential in colloidal suspensions [48j . 

Stochastic motions of a free particle moving along the 
streamlines of a sheared fluid and of a particle trapped 
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in the minimum of a potential, while exposed to a shear 
flow, have common characteristic signatures. Since the 
trapped particles are more suited for a thorough statis- 
tical analysis of its Brownian dynamics, we present cal- 
culations for particles trapped by a harmonic potential 
and exposed to either a linear shear flow or to a plane 
Poiseuille flow. Our analytical results show that shear 
flow causes characteristic signatures in the time depen- 
dence of the cross-correlation between particle displace- 
ments along orthogonal directions in the shear plane as 
well as an inclined elliptical particle distribution. Part of 
these results have already been applied and confirmed in 
a recent experimental study [TBj . 

For our calculations we utilize a Langevin model for 
the particle motion, where stochastic forces with differ- 
ent statistical properties may be used. Stochastic forces 
acting on suspended Brownian particles are caused by ve- 
locity fluctuations of the surrounding fluid. In a quiescent 
fluid the fluctuations of orthogonal velocity components 
are uncorrelated in the bulk [49j. Assuming such uncor- 
rected fluid- velocity fluctuations, and therefore uncorre- 
lated stochastic forces in the particle Langevin model, we 
show how shear flow induces cross-correlations between 
particle fluctuations along orthogonal directions. Con- 
versely, we show how the amplitudes of the stochastic 
forces acting on the particle can be determined by mea- 
suring the static correlations of the particle fluctuations. 

In section [II] the model equations of the Brownian par- 
ticle motion and their formal analytical solutions are pre- 
sented. The static correlation functions for the parti- 
cle's position and velocity fluctuations are derived in sec- 
tion where also the corresponding distributions are 
calculated in terms of the static correlations and under 
the assumption of Gaussian particle fluctuations. In ad- 
dition, the ratio between the principal axes of both dis- 
tributions are determined as well as the angle enclosed 
by each major axis and the flow direction. In Poiseuille 
flow the second derivative of the flow profile, as well as 
the fluctuations perpendicular to the flow lines, cause 
a shift of the particle's mean position in the potential 
via Au(r) ^ . The latter contribution is usually not 
taken into account, if the effective particle radius is de- 
termined via Faxen's law from the particle displacement. 
In section [IV] we present and discuss these results for 
the special cases of linear shear flow and plane Poiseuille 
flow. In addition we compare the analytical results with 
numerical simulations of the Langevin equation given in 
section [TT] and we suggest experiments to measure some 
of the flow-induced effects determined in this work. The 
article closes with a discussion and possible further ap- 
plications in section IVl 



II. EQUATIONS OF MOTION AND THEIR 
SOLUTIONS 

We consider a Brownian particle of mass m and effec- 
tive radius R suspended at the position r = (x, y, z) in 



a flow field with parallel streamlines in the x-direction, 
u(r) = u x (y)e x . We assume a velocity field 

u x (y) = (a + by + cy 2 ), (1) 

which corresponds for b = c = to a uniform flow, 
for a = c = to a linear shear flow with shear rate 6, 
and for c = — p-, b = to a plane Poiseuille flow be- 
tween two parallel walls at a distance 21. The particle 
is trapped by a harmonic potential with its minimum at 
r = (x ,yo,z ) = (0,0,0), 

U(r) = \, 2 . (2) 
The resulting linear restoring force is given by 

F p = -W = -/cr, (3) 

in terms of the force constant k. Such a potential act- 
ing on a colloidal particle may be realized by an optical 
tweezer [25| . 

A particle moving with the speed v = r in a flow of 
velocity u experiences, according to Stokes' law, a hy- 
drodynamic drag force F h = 67rr]R(u — v) proportional 
to the effective radius R, to the shear viscosity r\ and to 
the difference u — v between the velocity of the particle 
and the local flow velocity [49|, [50| . If the flow velocity 
u(r) is a nonlinear function of the spatial coordinates, as 
in the case of a plane Poiseuille flow in Eq. (pQ), one has 
according to Faxen's theorems [5l| an additional contri- 
bution to the drag force. This contribution includes the 
Laplacian of the velocity field and has in terms of the 
Stokes friction coefficient, ( = 6ttt]R, the form: 

F h = -c(v- u(r) - |Au(r)j . (4) 

The Laplacian of the flow field in Eq. (pQ), with the ab- 
breviation a = a + ^-c, gives the following expression for 
the hydrodynamic drag force, 

F h = -(r + ((a + by + cy 2 )e x . (5) 

The stochastic motion of a Brownian particle is caused 
by the fluctuations u of the fluid velocity-field u(r,t) 
[521 ]. The effects of u on a particle can be taken into 
account in a Langevin model by a random force F b (t). 
In uniform flows, namely with b = c = in Eq. (pQ), 
the cross-correlations of the velocity fluctuations of the 
fluid, (uiiij), lead to a vanishing cross-correlation of the 
random forces, (F b F b ) = 0, (i ^ j). The shear-induced 
contributions to U{ and F\ are a matter of current re- 
search [13, El. 

In our model we assume a Gaussian distribution of 
F b (t) with vanishing correlation time and mean value 

(Ft(t)Ff(t , ))=f ij 5(t-t) &ndi,j ex,y,z. (6) 
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For the moment we leave the fluctuation matrix 
unspecified, except to note that, according to time- 
translation- and time-reversal-invariance, it is symmet- 
ric. In uniform flows, namely with b = c = in Eq. (pQ), 
the matrix is diagonal with fa = 2ksT( [53| as men- 
tioned before, whereas in a shear flow the magnitude of 
the non-diagonal elements of depends on the shear 
rate, but shear- induced contributions are expected to be 
small [l3|, [H, [23j| . However, we allow non-diagonal ele- 
ments of fij for the moment to show in section IIVI how 
these non-diagonal elements may be determined by mea- 
surements of the velocity fluctuations of the particle. 
The net force acting on the particle, 



F /l + F^ + F 6 (t) : 



(7) 



together with Newton's law gives the Langevin equations 
of motion for the translational degrees of freedom of the 
particle 

mx = -(x -kx + ((a + by + cy 2 ) + F%(t) , (8a) 
my = -(y-ky + F*{t), (8b) 



rnz 



-(z - kz + F h z (t) . 



(8c) 



Introducing the vectors X = (x,v x ), Y = (y,v y ) and 
Z = (z,v z ) one may express the second order differential 
equations (|HJ) in terms of a system of coupled first order 
equations 



X = LX 



Y = LY 



Z = LZ 



rn 
m 



e v + 2/3(a + by + cy 2 )e v 



Herein we have introduced the matrix 

1 

J 1 -2/3 



(9a) 
(9b) 
(9c) 

(10) 



the damping constant f3 = 7^, the squared frequency 
uj 2 = ^ and the velocity unit- vector e v = (0, 1). The 
rotational motion of the particle, which may provide fur- 
ther corrections to the leading order fluctuation effects 
discussed here, is not taken into account in this work. 

The general solutions of the equations of motion (j9]) in 
terms of the initial conditions X(0), Y(0) and Z(0) are 
given by 



X(t) = e Lt X 



;(o) + f 

Jo 



Ht-r) F'Ar) 



dre ut -^ 



m 



(11a) 



+ 2/3/ dr e^-^ (a + by + cy 2 ) e v , 
Jo 

Y(t) = e Lt Y(0) + f dr e L ^ e v , (lib) 

Jo m 

Z(t) = e Lt Z(0) + / dr e^-^ ^1 e v . (11c) 
Jo m 



This system of coupled equations is the starting point of 
the following analysis, where the statistical properties of 
the particle's motion are characterized by the correlations 
of its position and velocity. 



III. DISTRIBUTION FUNCTIONS AND STATIC 
CORRELATIONS 

By taking the averages of Eqs. (fTTj) and using the van- 
ishing mean of the stochastic forces in Eq. (|6|) we see that 
the mean velocity of a particle in a harmonic potential 
vanishes, (v) =0. In the case 1*0 = the mean value of 
the deviations of the particle's position from the poten- 
tial minimum in the directions perpendicular to the flow 
vanish too, (y) = (z) = . However, the mean particle 
displacement in the direction of the flow given by Eq. (pQ) 
is non-zero: 



/ \ C - , fyy 
{x)= k a+ 2k^ C - 



(12) 



This equation is discussed in more detail in section HVl for 
specific flows. Since the spring constant k of the potential 
enters in different powers in equation (fT2|) . this formula 
might be employed for the experimental determination 
of the effective fluctuation magnitude f yy . 

With a combination of the coordinates of the parti- 
cle and its velocities to a single 6-dimensional vector, 
q = (r, v) = (x,y,z,v x ,v y ,v z ), the probability dis- 
tribution function of the particle, V(q), may be for- 
mulated in a compact form in terms of the deviations 
q = q — (q) = (x, y, z, v x ,v y ,v z ) from the mean value (q). 
If the fluctuations q are linear functions of Gaussian dis- 
tributed stochastic forces, which is the case for uniform 
and linear shear flows, they can be expected to be them- 
selves Gaussian variables [54| and can be described by a 
Gaussian distribution as follows: 



1 



a/^tt) 6 det(C) 



exp 



1 



cfC^q 



(13) 



Here the covariance matrix C = (qq T ) is used, which 
includes second order moments for the coordinates and 
the velocities at equal times (static correlations). The 
magnitudes of the elements dj depend on the correla- 
tion amplitudes of the stochastic forces, fy, and can be 
measured in experiments. Consequently, one may recon- 
struct the stochastic forces from the measurements, as 
discussed later. 

For Poiseuille flow the relations between particle fluc- 
tuations and stochastic forces are nonlinear due to the 
contribution oc cy 2 in Eq. (|lla|) . Therefore, the particle 
distribution V(q) is not necessarily Gaussian. However, 
in the course of this work we use the formula in Eq. (fT3|) 
also for particles in a potential, which are exposed to a 
Poiseuille flow, but with the covariance matrix C^- now 
determined in terms of the parameters for the Poiseuille 
flow. The validity of this heuristic approach will be tested 
later in Sec. HVBl bv numerical simulations. 
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A. Covariance matrix and angular momentum 



equations: 



The covariance matrix C consists of four 3x3 subma- 
trices 



C 



(14) 



that describe the autocorrelation for the positions C rr 
and velocities C vv and their cross-correlations C vr and 
C rv at equal times. The matrices 



p _ pi 




(15a) 
(15b) 
(15c) 



may be calculated in terms of the expressions given by 
Eqs. CP. 



Under the Gaussian assumption, the four-point corre- 
lations of the stochastic forces, which occur due to the 
quadratic contribution y 2 of the flow in Eq. (jllap , can be 
decomposed into two-point correlation functions accord- 
ing to Wick's theorem p54j : 



(F»( tl )FHt 2 )F b k (t 3 )F[>(n)) 



= fijfklS(t2 

+ fikfjiS(t 3 ■ 

+ filfjkS(t4 - 



■h)S(t 4 
h)S(t 3 - 



■t 3 ) 
ts) 

t 2 ). (16) 



Introducing the relaxation time of the particle in a har- 
monic potential 



c 

k ' 



(17) 



the non-zero components of the covariance matrices may 
be written in terms of the amplitudes of the stochastic 
forces and the flow parameters by the following set of 



(ft) 



1 fy 



4m 2 j3 ' 

1 fzz 

4m 2 /3' 



(z>) = 



(vl) 



I ~ ~ \ I ~ ~ \ 1 fxy 

( VxVy ) = (VyV X ) = -— ^ , 

m = (p) = + W P (f) 



(XVy) 



-(yv x ) 

i fx 



W 2 2 



- + -b (v ) H 

4m 2 ^ 2 W ; 3(l + 2r 2 w 2 ) 



(3 

LO 



+ T p b(xy} + -Ty(f) 2 . 



(18a) 
(18b) 
(18c) 

(18d) 
(18e) 

(y 2 ) 2 , (isf) 

(18g) 



The non-diagonal elements of C VVl namely the correla- 
tions of the velocity fluctuations of the particle, (v x v y ) 
and (v y v x ), are directly proportional to the amplitude f xy 
of the fluctuations given by Eq. (j6]). In contrast to this 
result for a fluctuating particle in a potential, one finds 
for a free particle in shear flow finite values for the cross- 
correlations (v x v y ) and (v y v x ) even in the case f xy = 
[l8| ]. However, as shown appendix lAl in the presence of 
a harmonic potential these correlations decay on a time 
scale 1/(2/3), which is very short for an overdamped par- 
ticle motion. In the case of a weak laser tweezer po- 
tential the particle relaxation time r p = Q/k becomes 
rather large and one obtains according to Eq. (jA8j) and 
Eq. (|A10j) for f xy = a contribution to the velocity 
correlation (v x (t)v y (t)) oc exp[— 2t/r p ] f yy , which de- 
cays slowly and gives in the limit of a vanishing poten- 
tial (k —> 0,r p — >> oo ) a constant contribution to the 
velocity- velocity correlation, which agrees with that in 
Refs. [18[ for a free particle, c.f. appendix lAl Therefore, 
finite values of (v x v y ) and (v y v x ) measured for particles 
trapped in a potential are a direct indication of cross- 
correlations of the stochastic forces along orthogonal di- 
rections, (F^F^) ^ 0. 

For 6/0 the contributions to the non-diagonal ele- 
ments (xy) and (yx) of the positional correlations C rr , 
which are related to f xy oc 6, are expected to be small 
HHH. For this reason f xy is neglected in the following. 
Both positional cross-correlations (xy) and (yx) are pro- 
portional to the local shear rate b and to the fluctuation 
strength f yy in the ^/-direction, cf. Eq. (|18d|) . 

The autocorrelation of the velocity fluctuations (vt\) 
in the flow direction includes several contributions. It 
depends on the shear rate, the second derivative of the 
flow and the fluctuation strength f xx and f yy . The cross- 
correlations between velocity and position appear only 
in the shear plane and they are proportional to the local 
shear rate. 

The sub-matrices C vr and C rv , describing the cross- 
correlations between the positional and the velocity fluc- 
tuations, are related to the mean angular momentum of 
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the particle: 

((rar x v) z ) m(xv y - yv x ) = -mb(y 2 ) . (19) 

B. Distribution of position and velocity 

Integrating the particle distribution function V(q) in 
Eq. (fT3|) with respect to its velocity degrees of freedom 
one obtains the particle's positional distribution function 
'P(r), which may be expressed in terms of the covariance 
matrix C rr as follows: 



P(i) 



1 



v /(27r) 3 det(C rr ) 



exp 



— -r C rr r 



(20) 




FIG. 1: The elliptical particle distribution in the xy shear 
plane is sketched for a non-zero shear rate, di and di are the 
two principal axes of an ellipse along which V(y) is constant. 
The Gaussian profiles along the two principal axes are also 
indicated. 

A sketch of the distribution ^(f) is shown in Fig. [T] for 
a linear shear flow b ^ and a = c = 0. The non-diagonal 
elements of the symmetric matrix C rr , cf. Eq. (|15aj) and 
Eq. ([T8]h describe the cross-correlations of the particle 
fluctuations in the x- and ^/-direction and therefore the 
deviation of ^(f) from a spherically symmetric distribu- 
tion to an ellipsoidal one in the shear plane. The principal 
axes of the particle's positional distribution are given by 
the eigenvectors w p; i ? 2,3 of the corresponding eigenvalues 
Cp ; i,2,3 of the matrix C rr . 

In the z-direction one has c p -s = (z 2 ) and w p; 3 = 
(0,0,1). For the other two directions in the xy shear 
plane one obtains 



Cp;l,2 = \ ((X 2 ) + (f)) 



(21) 



Wp;l 



±\\jm) 2 + ({* 2 )-{y 2 ))\ 



(xy) 



\(xy)\ V ' (xy) 



• 



With the orthogonal transformation matrix 

D = (Wp;l,W p;2 , W p;3 ), 



(22) 



the f-dependence of the exponential function ()20|) can be 
rewritten as follows 



V c- 



l - (D T f) T diag(c-i,c-i,c-i) (D T f) . 

(23) 



The eigenvalues c p -i and c p; 2 determine the length of the 
principal axes d P] i^ = yjc^i^ along the directions w p; i 
and w p; 2 of an ellipse, where the longer one is rotated 
counterclockwise with respect to the x-axis by an angle 
p , which is given by the expression 




(24) 



(x 2 ) - (y 2 ) 

(xy) 



The ratio between the minor and the major axis is given 

by 



Vn 



Cp;2 
Cp;l 



(25) 



(x 2 ) + (y 2 ) - V4(iy) 2 + ((x 2 ) - (y 2 )) 2 



V (x 2 ) + (y 2 ) + vW+PFW / 

Analogous to the particle distribution given above, one 
can also derive an expression for the probability distri- 
bution V(y) of the particle's velocity, which is obtained 
by integrating out the positional degrees of freedom in 
Eq. m 



V{v) = J drV(r,v) 



V(2^) 3 det(C„) 



exp 



(26) 



The eigenvalues and eigenvectors w u; 1,2,3 of C vv 

are determined in a similar manner as for C rr . Again 
the eigenvalue c v; s = (v 2 ) and the principal axis w v; s = 
(0,0,1) perpendicular to the shear plane are obvious. 
The remaining two eigenvalues and eigenvectors are: 



C-u-1,2 



w v; i 



\ m + (vD) 



±yMv x v y ) 2 + ((v 2 )-(v 2 y )) 2 , (27) 




(28) 
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They have the same structure like those for the covari- 
ance matrix C rr . However, while the non-diagonal ele- 
ments of C rr are directly proportional to the local shear 
rate 6, the non-diagonal elements of C vv and therefore 
the angle enclosed between the principal axis of the dis- 
tribution of the velocity fluctuations and the x-axis, 



tan 



(VxVy) 




(29) 



(v 2 y y 



(VxVy) 



depends only weakly on the local shear rate, since the 
contribution of b to the shear-induced cross-correlation 



f xy is small in magnitude, compared to fa [14],[23j. From 
a measurement of the distribution of the particle's veloc- 
ity one may calculate the eigenvalues c v; i^ and the angle 
<p v , which may enable the determination of the stochastic 
force correlations f xx , f yy and f xy . 

The mean kinetic energy of a trapped particle in shear 
flow is composed of two contributions, one induced by 
the distribution of the fluctuations and an additional one 
by virtue of the rotational part of the flow: 



L , m 2 



1 fxx H~ fyy H~ f , 

4 m 2 P 



1 + 2r p 2 w 2 



(30) 

(f) 2 ■ 



If the correlations of the stochastic forces are assumed to 
be independent of the flow, the mean energy of the par- 
ticle will be increased by the flow, since all terms in this 
equation are positive. But without knowing the explicit 
expressions for the force correlations, it is not clear how 
the energy is really changed. Recent calculations show 
that the mean kinetic energy of a fluid without an im- 
mersed particle increases in a shear flow [55, 56]. This in- 
dicates that the stochastic forces and the particle's mean 
energy, Eq. ([30]) . may be amplified as well. 



of shear rate 7, we may insert in Eq. (pQ) y = yo + y, 
where y describes the deviation from y$. Consequently, 
in terms of the coordinates y, the three coefficients in 
Eq. (pQ) take the following form 



6 = 7, 



a = a = 72/0 



and c = . 



(31) 



In the case 1/0 / the flow velocity has a finite value 
at the center of the trap resulting in a non-zero mean 
position of the particle in flow direction, which is given 
via Eq. (fT2|) by the formula 



(x) = jr p y . 



(32) 



With the identifications (f3Tj) the elements of the co- 
variance matrix C for the particle fluctuations around 
the potential minimum may be determined for a linear 
shear flow by Eqs. ([T8|) in terms of the noise amplitudes 

fxxi fxy and fyy 

An interesting question is, how to detect the noise am- 
plitudes fij in terms of the measured auto- and cross- 
correlations of the particle's position- and velocity fluc- 
tuations. According to Eq. (|18cj) there is a direct re- 
lation between the velocity correlation (v x v y ) and the 
noise magnitude f xy . Therefore, a direct measurement 
of this velocity correlation, if experimentally possible, 
would give f xy , or, the other way around, a non-zero 
value of f xy ^ is required in order to obtain non-zero 
values of the cross-correlations between these orthogonal 
velocity components. 

As already mentioned in the introduction, it is favor- 
able to investigate the particle fluctuations around a po- 
tential minimum rather than those of free particles since 
trapped particles can be investigated over a long period 
of time, as demonstrated by several experiments, see e. 
g. Ref. [36j. This opens the opportunity for the determi- 
nation of the magnitudes of the noise fij in terms of the 
positional fluctuations via Eqs. (fl8|) . which in the case of 
a linear shear flow take the following explicit form 



IV. RESULTS FOR SPECIFIC FLOWS 



fyy = 2k((f) , 

f xy = 2k( ( (xy) - \ir p {y 2 ) j , 



(33a) 
(33b) 



The results concerning the effects of a linear shear and 
a plane Poiseuille flow on the Brownian motion of a parti- 
cle in a harmonic potential are presented in this section. 
Namely the properties of the correlation functions and 
the particle's distribution are analyzed as a function of 
the flow parameters. They share common trends for both 
flows but there are also some characteristic differences, 
which are described in this section. 



A. Linear shear flow 

If the center of the harmonic potential at ro = (0, 7/0 > 0) 
does not coincide with the center of the linear shear flow 



2kcl(x 2 )-^(f)-jr p (xy) 



(33c) 



In the case of isotropic stochastic forces with negligible 
cross-correlations, f xx = f yy and f xy oc (v x v y ) « 0, one 
has in a linear shear flow still a non- vanishing positional 
cross-correlation (xy) = / yr p (y 2 )/2. Its magnitude is de- 
termined by the shear rate 7 and the noise strength via 
(y 2 ). According to this behavior and because (x 2 ) ^ (y 2 ), 
we expect an anisotropic distribution of the positional 
fluctuations V(r), as sketched in Fig. [T] and as discussed 
below. For f xy ^ the anisotropy of the positional dis- 
tribution has an additional contribution that depends on 
the magnitude of f xy . 
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During the rest of the present section IIV Al inertia 
effects are neglected and in addition, we assume an 
isotropic noise distribution with f xy = f yx = and 
fxx = fyy = fzz = 2k B T(. Both are good approxi- 
mations for many experiments focusing on leading order 
effects of a shear flow on the fluctuations of particles. 
Taking into account the time-dependence of the posi- 
tional fluctuations, as given by Eqs. ([TT]) . one obtains 
in terms of the dimensionless Weissenberg number, 



Wi = r p j , 



(34) 



and the Heaviside step function Q(t) the following time- 
dependent correlations: 



(35a) 
(35b) 



w)m) = (mm) = ^e- w/r ', 

(mm) = ™ (l + 2^ 8(f)) e-l* 1 ^ . (35c) 



are in agreement with experimental data as described re- 
cently in Ref. jl5|. The expression (x(t)y(0)) is propor- 
tional to the Weissenberg number and takes its maximum 
at half of the particle's relaxation time tm ax = r p /2. For 
the correlations of the velocity fluctuations of the fluid in 
orthogonal directions a similar signature as in Eq. (|35c|) 
has been found [l3|, where however the mechanism is 
slightly different. 

A comparison of the absolute values of the particle's 
fluctuations with experimental results may be difficult. 
However, one obtains from Eqs. ([35]) and in terms of the 
dimensionless Weissenberg number the following normal- 
ized ratios of the static correlations 



(mm) = wi 

2 



(mm) 
(mm) 

(5(0)5(0)) l + Wi 2 /2 



Wi/2 



(36a) 
(36b) 



The left hand side and the right hand side of Eqs. (f36|) 
can be measured independently in different experiments 
and the results can be compared afterwards. 



The cross-correlation between fluctuations along orthog- 
onal directions, as given by the last equation, is shear- 
induced and its asymmetry (x(t)y(0)) ^ (x(0)y(t)) with 
respect to time reflection t — »• — t is one of the important 
effects of shear flow on the distribution of fluctuations. 




2 

V 
-2 
-4 
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FIG. 3: Positional distribution of a particle in a harmonic 
potential with its minimum at 1*0 = (0,0,0) and exposed to a 
linear shear flow, as obtained by a stochastic simulation of the 
Eqs. JS} in the over-damped limit for the Weissenberg number 
Wi = 2 and the relaxation time r v = 1. 



FIG. 2: Time dependence of the shear-induced cross- 
correlations (x(t)y(0)) (dashed line) and (x(0)y(t)) (solid line) 
as given by Eq. (|35cp for the Weissenberg number Wi = 1 and 
the relaxation time r p — 1. 

For t > 0, the algebraic prefactor in Eq. (|35cp illus- 
trates that a fluctuation y(0) ^ of a particle is carried 
away by the flow along the x-direction before the initial 
displacement starts to relax remarkably. This leads, dur- 
ing an initial period of time shorter than the relaxation 
time r p , to the growth of (x(t)y(0)), while (x(0)y(t)) 
decays monotonically (solid line in Fig. [2] correspond- 
ing to t < 0). The predicted elementary signatures for 
the shear-induced cross-correlations, as shown in Fig. [21 



An anisotropic distribution of the particle's velocity 
P(v), as given by Eq. ([26|) . is only obtained in the case 
of a finite cross-correlation f xy of the stochastic forces. 
In contrast to this, the particle's distribution V(r) has 
in the overdamped limit an elliptical shape in the xy 
plane, even for a vanishing cross-correlation magnitude 
fxy = 0. In this limit an elliptical distribution V(r) is 
shown in Fig. [3] for the Weissenberg number Wi = 2 and 
the relaxation time r p = 1. In Fig. [3] the time evolution 
of the particle's position, obtained by simulations of the 
basic equation (|HJ), is plotted at equidistant times. This 
positional distribution can be characterized by the an- 
gle <p p of the ellipsoid and the ratio V p of its principal 
axis. The general expressions for <p p and V p , as given by 
Eqs. (|2l|) and ([25]) . can be further simplified in the case 
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of a linear shear flow to functions of the dimensionless 
Weissenberg number Wi only: 



tan < 




(37a) 
(37b) 



The two expressions in Eqs. (|37j) suggest measure- 
ments of the shear- induced particle fluctuation effects, 
which are complementary to the measured static corre- 
lations. In experiments the particle positions may be 
recorded at equidistant times. By plotting these sub- 
sequent particle positions in the shear plane, a similar 
distribution is expected as shown by our numerical simu- 
lation in Fig. [3l From such an experimentally measured 
distribution for different shear rates the angle (j) p and the 
ratio between the principal axes, V p , may be determined. 
If the determination of the Weissenberg number Wi is 
difficult or if the precision is not sufficient Eqs. (f36|) and 
(|37|) allow a consistency check between different aspects 
of the particle fluctuations, without a separate measure- 
ment of Wi. A cross-check has recently been performed 
in an experiment in which a good agreement between 
both approaches has been found, cf. Ref. [l5| . 




V p 0.6 




2.0 2.5 



FIG. 5: The ratio V p of the principal axis of the ellipse in 
Eq. (|37b[) , as a function of the Weissenberg number Wi. The 
triangles are the result of a simulation of Eq. JS} in the over- 
damped limit. 



The triangles in those figures have been obtained by nu- 
merical simulations of Eq. ((£)) for the same values of the 
Weissenberg number as used for the analytical curves. In 
our simulations of the Langevin model we have assumed 
an isotropic and Gaussian distributed white noise and 
vanishing cross-correlations f xy = 0. The simulation re- 
sults confirm the assumption of a Gaussian distribution 
of the particle positions in the case of a shear flow in 
Sec. IIII1 The analytical and numerical results on (j) p and 
V p in Fig. [Hand in Fig. [5] are in good agreement. 



B. Plane Poiseuille flow 

A Brownian particle is trapped by a harmonic potential 
close to i*o = (0,?/o,0) and exposed to a plane Poiseuille 
flow along the x-direction, 



u (y) 



u, 



(38) 



FIG. 4: The angle of the ellipse in the xy plane with respect 
to the flow lines as a function of the Weissenberg number Wi, 
as it is given by Eq. (37aj). The triangles are the result of a 
simulation of Eq. JSJ) in the over-damped limit. 

In the limit of a vanishing Weissenberg number, Wi — >• 
0, the angle (j) p of the ellipsoidal particle's distribution 
tends to <p p = 7r/4 and the positional variance becomes 
isotropic, (x 2 ) = (y 2 ), corresponding to the ratio V p = 1. 
This trend is similar to the dependence of the orientation 
of vesicles in shear flow [Hzl or to the local orientation of 
the order parameter of nematic liquid crystals in plane 
shear flows (58) . 

If 7 and therefore Wi is increased, the angle (j) p be- 
tween the longer semi- axis and the x-axis as well as the 
ratio V p will decrease as shown in Fig. [H and in Fig. [5j 



The particle fluctuations x and y describe deviations with 
respect to the mean values yo and (x), which is not zero 
as determined below. The flow velocity along the re- 
direction may be expressed for further analysis in terms 
of the fluctuations y as follows 



1 _ yl 
i 2 



u p y „ Up „ 2 



I 2 



■V-pV 



(39) 



Comparing this expression with Eq. (pQ) the coefficients 
in the latter equation are given by 



a ■ 
b 



>p ~ XVo ■ 



(40) 
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Here x describes the second derivative of the velocity 
profile, b the local shear rate and 



Eqs. ([15]) . then reduce to 



Wi(2/ ) := -2r P x^/o, 



(41) 



the local Weissenberg number Wi(yo)- With these iden- 
tifications the elements of the covariance matrix C are 
again given via the expressions in Eqs. (fTHj) in terms of 
the strength of the noise, f xx , f xyi f yy and f zz and the 
flow parameters. 

The mean position of the particle in a plane Poiseuille 
flow can be determined via Eq. QX 



(x) = T p 



X 



x(y 2 ) 



(42) 



In contrast to a linear shear flow, it includes a contri- 
bution depending on the particle's radius R, which is a 
pure deterministic effect due to Faxen's theorem [2I. lElj. 
The last term on the right hand side describes an addi- 
tional shift based on the positional variance (y 2 ) in the 
^-direction. Both contributions are proportional to the 
second derivative \ °f the flow profile and are therefore 
not present in linear shear flows. 

In experiments Eq. ([42]) may be used to measure the 
spatial variation of the flow profile by detecting the mean 
displacement of a particle of radius R out of the opti- 
cal tweezer potential. For such a measurement usually 
deterministic formulas are used to describe the relation 
between the displacement and the flow velocity. But 
Eq. ([42]) indicates that a correction due to thermal mo- 
tion has to be taken into account. 

The relations given by Eqs. (fT5]) relate the fluctuations 
of the velocity and the position of a particle to the ex- 
ternally controlled flow properties and the magnitudes of 
the thermal fluctuations. Consequently they allow, in a 
similar manner as to the linear shear flow, a determina- 
tion of the magnitudes of the stochastic forces in terms 
of the measured covariances of the particle fluctuations: 

(43a) 

;Wi(y ) if)) , (43b) 



2 T p UJ Z 

2r^ 



fyy 


= 2k({y 2 ) , 


fxy 


= 2K({xy) 


f XX 


= 2kC^(x 2 ) 








3 P x 1 1 



(y 2 y-Wi(y ){xy)j . (43c) 

The difference compared to Eqs. ([33]) is an additional 
contribution to f xx , which depends on the second deriva- 
tive of the Poiseuille flow. Further comments made above 
for a linear shear flow hold as well. 

Similar to the end of the previous section IIV A[ we 
neglect from here on the particle inertia and we assume 
in addition isotropically distributed noise with f xy = 
and fa = 2ksT( > . The static correlations of the parti- 
cle's positional fluctuations in the shear plane, given by 



(x 2 ) 



if) 



k B T 
k 

k B T 



(xy) = (yx) = 



:Wi(7/ ) 2 



k B TWi(y ) 
k 2 



2 2 k BT 



(44a) 
(44b) 
(44c) 



The static cross-correlation (xy) has the same depen- 
dence on the Weissenberg number as in the linear shear 
case Eq. ()35cj) at t = 0, if the local shear rate of the 
Poiseuille flow at the potential center is taken. In con- 
trast to a linear shear flow profile, the mean-square 
displacement of the particle in x-direction, given by 
Eq. (j44ap , includes an additional constant contribution, 
G := 2/3 (xTpksT '/k) 2 , which is independent of the posi- 
tion of the potential minimum in the Poiseuille flow and 
therefore also independent of the local Weissenberg num- 
ber. Besides the dependence on the second derivative x 
of the flow profile, the contribution G is a direct func- 
tion of the fluid temperature. Due to this contribution 
in Eq. ()44aj) one has (x 2 ) ^ (y 2 ) and the isotropy of the 
particle's positional distribution in the shear plane is bro- 
ken. This broken rotational symmetry also changes the 
analytically determined distribution function in Sec. lIIIBl 
along with the correlation matrix for a Poiseuille flow. 

In addition to the calculations, we have performed sim- 
ulations of the particle dynamics where we used isotropic 
and Gaussian distributed white noise in the overdamped 
version of Eq. (|8|). The resulting distributions for the 
particle's position are shown in Fig. [6] for three different 
positions yo of the potential center. With the potential 
minimum at the center of the flow (yo = 0) the numerical 
results show a broken mirror symmetry in x-direction of 
the particle's positional distribution, c.f. Fig. JHfci). One 
recognizes a parachute shape that is similar to the well 
known conformation of vesicles and red blood cells in 
the center of a Poiseuille flow [59|. If one now makes 
the heuristic assumption of a Gaussian distribution for 
the particle's position in the case of a Poiseuille flow, 
as in Sec. (|III|) . but with the correlation matrix dj de- 
termined in terms of the Poiseuille-flow parameters, one 
expects for the parameters in Fig. [6^) an elliptical shape 
of the particle distribution. Indeed, the ratio between 
the principal axis in Fig. [6^) is slightly smaller than 1.0. 
But within this analytical approximation, the ±x sym- 
metry is not broken, which indicates the limitation of the 
heuristic approach. 

Away from the center of the Poiseuille flow, for finite 
values of yo ^ 0, the ±y symmetry of the particle's po- 
sitional distribution is also broken, as shown in Fig. [Bb) 
and in Fig. [6fc). With increasing values of yo the lo- 
cal shear rate acting on the particle increases as well 
as the local Weissenberg number Wi(?/o)- Consequently 
the cross-correlation (xy) in Eq. (|44cj) becomes non-zero 
and the particle's positional distribution in the xy plane 
approaches, according to our analytical results, an el- 
liptic shape as indicated by the ellipses in Fig. Eb) and 




FIG. 6: Positional distribution of a Brownian particle in the xy plane captured by a harmonic potential with its minimum at 
xo = zo = and for different values of yo. The particle is exposed to a Poiseuille flow and its dynamics have been calculated 
by stochastic simulations of the Eqs. ([8} in the over-damped limit for x = 1, ksT = 0.1 and r p = 1. 



Fig. [6t). Again the full numerical simulations show de- 
viations from the elliptical shape. 

However, the inclination of the distribution and the 
inclination of the analytically determined ellipses agree 
rather well and therefore a determination of the angle 
(j) p for a Poiseuille flow according to Eq. ([24]) , similar to 
that in the previous chapter, is reasonable and c/) p has the 
following form 



at yo-,max, as given by the following equations: 



tan ( 




Wifoo) 



Wi(yo) + - 




3Wi(2/ ) 



We have shown in Sec. IIV A\ how the angle <p p and the 
ratio V p depend on the Weissenberg number of the linear 
shear flow. Since the local shear rate in a Poiseuille flow 
depends on the location yo of the minimum of the po- 
tential, one may plot (j) p and V p (which is calculated via 
Eq. ([25]) ) as a function of yo as shown in Fig. Since 
the particle's distribution is anisotropic in the xy plane 
even for y = 0, see Fig.[6fci), the corresponding principal 
axes are always unequal. This behavior is reflected in 
Fig. 03) where the inequality V p < 1 holds for all values 
of yo- The dependence of the angle <p p on the local shear 
rate in a Poiseuille flow differs from the case of a linear 
shear flow: the function (j) p (yo) is always well defined, 
even at the center of the flow, where it vanishes. This is 
one consequence of the asymmetry of the particle distri- 
bution. With increasing values of \yo\ and therefore with 
increasing values of the local shear rate, the angle \(j) p \ 
increases as well until it reaches some maximum value 



yo Vl 



± 



3 



=F arctan 




r2 x 2 (y 2 ) - 2t p x 




(46a) 



(46b) 



yo,max depends only on (y 2 ) and therefore on the width 
of the particle's distribution in ^/-direction, which is de- 
termined by the ratio between the thermal energy and 
the spring constant related to the harmonic potential 
acting on the particle. Increasing \yo\ beyond \yo,max\ 
the local shear rate, Wi(yo), increases too, but (j) p starts 
to decrease; a similar behavior as seen in section IIV Al 
In the range \yo\ > \yo,max\ the local shear dominates 
the curvature effects more and more and the Poiseuille 
flow resembles a linear shear flow. In addition, the par- 
ticle's distribution approaches an ellipse as obtained by 
the heuristic approximation. As mentioned above, the 
heuristic analytical approach becomes exact in the case 
of a linear shear flow. 

In Fig. [7] the analytical results on (j) p and V p have been 
compared with our numerical simulations. In spite of 
the fact that we assumed for our analytical calculations a 
Gaussian distribution of the positional fluctuations x and 
y, the results show surprisingly good agreement. The ma- 
jor reason for this good agreement is that we calculated 
the inertia tensor for the particle distribution, which in- 
cludes only second order moments as assumed for the 
Gaussian distribution. 
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Besides the factor 2 in front of R 2 this expression corre- 
sponds to that in Eq. 



1 


1 


1 


i 



Vn 




FIG. 7: The angle cj) p in the xy plane according to Eq. (|35|) 
is shown in a) and the ratio V p between the two princi- 
pal axes , as determined by the Eqs. and Eq. ([25]) . is 
shown in b). Both are plotted as a function of the posi- 
tion of the potential minimum yo and for the parameters 
X = t p = 1, (y 2 ) = ksT/k = 0.1. The triangles are the 
result of simulations of the Eqs. (0 with the same parameter 
set. 



C. Pipe Flow 

For a flow profile in a pipe of radius <i, cylindrical polar- 
coordinates are appropriate. Similar to a plane Poiseuille 
flow the dependence of the velocity profile on the radial 
position p is also quadratic. Consequently, the p depen- 
dence corresponds now to the y dependence in the case 
of the plane Poiseuille flow. The velocity profile for the 
pipe flow, where the pipe axis coincides with the x axis, 
is of the following form, 



1 d 2 



(47) 



The results described above in section HV Bl apply quali- 
tatively also for the pipe flow. 

With the potential minimum located at po the coeffi- 

2 

cients in the Eq. (pQ) are a = u p (l — ^), c = —\ '•— — 
and b = 2poc. Moreover, one has to consider a = a+|P 2 c 
in Eq. (j5j). The major difference to the plane Poiseuille 
flow is the prefactor | instead of | . This is a consequence 
of the 2D Laplacian instead of a ID second derivative in 
Eq.(j4]). The mean-position of the particle is in terms of 
the pipe flow of the following form 



V. CONCLUSION 

In this work we have calculated analytically and nu- 
merically the autocorrelations and cross-correlations be- 
tween different components of the velocity and the posi- 
tional fluctuations of a Brownian particle in a harmonic 
potential, which is exposed to different shear flows. In ad- 
dition, the particle's probability distribution in the har- 
monic potential has been determined as a function of the 
flow parameters. By solving an appropriate Langevin 
model, cross-correlations between velocity and positional 
fluctuations along orthogonal directions have been found 
and several suggestions for experimental measurements 
are made. 

Cross-correlations between orthogonal velocity compo- 
nents occur only if there is already a cross-correlation 
between stochastic forces along orthogonal directions in 
the related Langevin model. On the other hand, we find 
cross-correlations between particle fluctuations along or- 
thogonal directions without cross-correlations of orthog- 
onal stochastic force components and their magnitude 
increases with the dimensionless shear rate, the Weis- 
senberg number. 

There are recent calculations on shear-induced cross- 
correlations between orthogonal fluctuations of freely 
floating particles [HI, [HI . They have the same origin as 
those discussed in this work. However, while measure- 
ments of the shear-induced cross-correlations of freely 
moving particles may be difficult, shear-induced cross- 
correlations for particles in a potential can be measured 
in a controlled manner. This insight was the basis of 
recent successful measurements on the cross-correlations 
between particle fluctuations along orthogonal directions 
as described in Ref. [T^j. In this experiment, harmonic 
potentials for small latex spheres are induced by opti- 
cal tweezers. The particle is simultaneously exposed to 
a linear shear flow in a special flow cell and its Brown- 
ian motion is investigated directly. In the same experi- 
ment and in a forthcoming work, the cross-correlations 
and anti-correlations between two particles captured in 
two neighboring potentials and exposed to shear flows 
are also investigated. 

The measurements on the pos itional cross-correlation, 
(x(t)y(0)), presented in Ref. [15j, exhibit a similar max- 



(x) 



XT P (y 2 



(48) 



imum as predicted by the expression in Eq. (|35cj) and as 
shown in Fig. [2l This maximum is a typical signature 
of Brownian motion in shear flow and it is found where 
the shear rate approximately equals half of the particle's 
relaxation time. Also the elliptic shape of the particle's 
distribution, P(r), as shown in Fig. 03 has been measured 
in Ref. [15j. The angle <p pi enclosed by the major axis of 
the distribution P(r), and the ratio V p between the length 
of the two principal axes on the one hand and the magni- 
tude of the static correlations on the other hand focus on 
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different aspects of the dynamics of a particle in a poten- 
tial exposed to a shear flow. The interrelation between 
these different aspects of the Brownian particle dynam- 
ics can be used for consistency checks in experiments; 
such cross-checks are described Ref. [15| for the case of a 
linear shear flow. They are also possible in experiments 
with Poiseuille flow as forthcoming measurements may 
demonstrate. 

Shear-induced correlations of fluctuations of the fluid 
velocity with respect to a linear shear flow were found 
theoretically [13, Hi. The y 

are traced back to the non- 
normal property of the Navier-Stokes equation linearized 
around the linear shear profile [l3| and these velocity 
fluctuations play an important role in a shear flow for 
its instability and the onset of turbulence. The cross- 
correlations between these velocity fluctuations exhibit a 
similar extremum as given by Eq. ()35cj) , but it is a slightly 
different mechanism leading to this similar behavior on 
the level of the correlation function. This interrelation 
may be discussed in more detail in future work. 

Stochastic forces acting on a particle in a fluid are 
caused by the velocity fluctuations of the fluid surround- 
ing a particle. In quiescent fluids the correlations for ve- 
locity fluctuations are isotropically distributed. For rea- 
sons of simplicity, this is quite often assumed in models 
to investigate the Brownian motion of particles in shear 
flow. In that case, cross-correlations between stochastic 
forces along orthogonal directions vanish in the Langevin 
equation of motion. Shear induced cross-correlations be- 
tween orthogonal components of fluctuations of the fluid 
velocity are indeed small compared to the shear inde- 
pendent contributions [3, [23| and hence this simplifica- 
tion is reasonable. However, to which extent such cross- 
correlations of the fluid velocity may quantitatively mod- 
ify the presented results on the cross-correlations of par- 
ticle displacements needs further investigations. 

In numerical simulations of the dynamics of a Brown- 
ian particle in a potential and exposed to a Poiseuille flow 
we found higher order correlations for the particle fluctu- 
ations of non- Gaussian behavior. However, for deriving 
our analytical results on the angle p , which the major 
axis of the particle distribution encloses with the flow di- 
rection, and the ratio V p between the lengths of the two 
principal axes, we assumed Gaussian distributed particle 
fluctuations also in Poiseuille flow. Accordingly, a perfect 
agreement between the results of the numerical simula- 
tions and the analytical calculations on (j) p and V p could 
not be expected. Nevertheless, we find good agreement 
between both approaches, cf. Fig. [7J especially further 
away from the center of the Poiseuille flow, where the lin- 
ear contribution in the flow profile dominates and where 
the assumptions are better fulfilled. 
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Appendix A: Cross-correlation (v x (t)v y (t)) 

In this appendix we discuss the correlation (v x (t)v y (t)} 
and we show that it includes, besides the contribution 
oc f xy in Eq. (j!8cj) . in the limit k — >• an additional 
contribution oc f V y , which corresponds to the result de- 
scribed in Ref. [18j. From Eqs. (JTTJ) one obtains 



Vy(t) 



dra vv (t — t) 



m 



2f3b y(r) 



dra vv (t - t) — , 

rn 

t ^ ( r ) 

dr a rv (r - r) - 



m 



with the abbreviations 



= el e L « 



a rv (t) = e L( ^ e v with e r 



,(A1) 
(A2) 
(A3) 

(A4) 
(A5) 



The equal-time velocity correlation takes then the form 

(V X (t)Vy(t)) = 



dr dr'a v 
Jo Jo 



,(t - r)a vv (t - r )( 



2f3b / dra vv (t - r){y{r)v y (t)) . 
Jo 



(A6) 



With the correlation between the position y(t) and the 
velocity v y (i), 



(y(r)vy(t)) 



ll 

rn 



mm(t,r) 

m dT"a rv (T-T")a vv (t-T"). 





(A7) 



one ends up with the expression 



2f3b^r [ dr [ df a vv (t - r)a rv (r - f)a vv (t - f) 
Jo Jo 

fxy Gi(t) 



2f3b 



fyy G 2 (t) 

m 2 4 



(A8) 



where we have introduced the functions G\(t) and G2(t) 
as well as the parameter 5: 



d(t) 



G 2 {t) = 



P cos 2 (St) cos(5t) sin ((ft) 



S 2 

-2/3 2 



25 2 p 



-2/8t 



s 
1 



(A9) 



sin(cft) cos(ft) /?sin 2 (<ft) £ sin 2 ((ft) 



2<5 3 S 4 
fit cos(ft) sin(ft) t cos 2 (St) 



25 2 



S=V~u 



2d 2 

e~^\ (A10) 
(AH) 
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In the limit t — »• oo the contribution G±(t) remains fi- 
nite, limt^ooGi(t) = jp. For finite values of r p = (/k 
the function G?2(t) vanishes on a timescale £ ^> 1/(2/3) 
too. But in the limit of a vanishing laser tweezer poten- 
tial, k ~ uj 2 0, one has 5 ~ ±i(/3 - cj 2 /( 2 /^) + •••) 
and in this limit some contributions of G2(t) cancel 
each other and the remaining terms are proportional to 
exp(— 2t/r p ). These contributions do not decay in the 
limit t p — >• oo and one gets G^t) = l/(8/3 3 ). In this case 
one obtains even for f xy = a non-zero cross-correlation 



(V X (t)Vy(t)) 

described in [!£ 



4C 2 



-, which is exactly the same result as 



Appendix B: Anisotropic trap 



In experiments the force exerted by the optical tweezer 
on a colloidal particle may be anisotropic in the shear 
plane. For instance, the force constant may be differ- 
ent along and perpendicular to the laser beam. If the 
laser beam does not hit the shear plane perpendicularly 
the difference between the force constants in x- and y- 
direction increases with the deviation from the orthogo- 
nal direction to the shear plane. Let's assume in Eq. d3j 
a different force constant for each direction: k x , k y , k z . 

The relaxation of the particle is now different along 
different directions around the potential minimum. In 
order to take this effect into account, we introduce three 
relaxation times in the same way as in Sec. IIII[ 



(Bl) 



and we further define the fraction 



/i := 



k x 



(B2) 



between the eigenfrequencies in the shear plane, which is 
used to express k y by k x via k y = fik x . 

The covariance matrix C from Sec. [TIT] depends now in 
a complex manner on the different force constants, but in 
the limit /i — >• 1 and k z = k x the result from the previous 
sections are obtained again. 

The different static correlations as given for the 



isotropic forces in Sec. HTT1 change as follows: 
<^> = 7lfe, (B3) 



(V 



1 fyy _ fty) 



4 m 2 j3uOy UJy ' 

1 ~ /i 2 2(V X Vy) 



\ fxx 1t2/~2\o i ^ 2 ^x T y;x 
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(B4) 

(^) 2 e 2 
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[ ' Am 2 pool to 2 ' 
where the following abbreviations have been introduced: 

,2 



(B6) 
(B7) 
(B8) 
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e 2 = 



t 2 . x u> 2 x + ±(1 + 2h 2 ) 
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Since the y- and ^-displacements are independent of 
any other direction, as can be seen in Eqs. (fTTj) . their 
autocorrelations are functions of the individual force con- 
stant only and do not depend on ji. However, the au- 
tocorrelations for the velocity and the position in in- 
direction depends in a complex manner on the differ- 
ent force constants. The same applies for the cross- 
correlations in the shear plane as follows 



(VxVy) 



1 fxy 
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